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Abstract 

We investigate a system of two weakly coupled oscillators. It is 
shown that an external periodic perturbation can lead to the capture 
into resonance. Asymptotic description and numerical simulations are 
presented. We have obtained the explicit formula for the threshold 
value for the amplitude of the perturbation that leads to the autores- 
onance. 



Introduction 

A system of nonlinear oscillators is often used as a standard model for oscil- 
lating processes. The oscillators have eigen-frequencies and can resonantly 
interact between each other or with an external perturbation. A slow varia- 
tion of frequencies can lead to the increase of amplitudes of oscillations. This 
phenomenon is usually called the capture into autoresonance and was found 
in PQ, [2], [3]. Later similar effects were observed in various fields of modern 
physics 0], 0,[6],[7]. 

Various aspects of autoresonance solutions were investigated for one- 
dimensional problems j8], [9] and for mult i- dimensional problems [20], |16j . 
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In this article we investigate the system of two weakly coupled nonlinear 
oscillators under a small perturbation 

x" + uPx = eaixy + e^expjic/?} + c.c), (1) 
y" + (2cu) 2 x = ea 2 x 2 , 

here ip = (u + ear)9, r = e9, e is a small positive parameter. 

The constants in the system have the following sense u> = const is fre- 
quency of the oscillator with the amplitude x, ct\ and a 2 are parameters of 
the nonlinear coupling, 7 is the amplitude of the external perturbation, a is 
the derivative of detuning of frequency of external perturbation with respect 
to slow time r. 

When e = the eigen- frequencies relate to each other as 1 : 2. If e 7^ 
then this relation leads to the parametric resonance. When external pertur- 
bation 7 = the nonlinear oscillators exchange energy [10] , [11] and the order 
of solutions does not change. 

Our goal is to obtain a solution with an increasing amplitude due to a 
small oscillating perturbation when 7 7^ 0. There is a standard way to obtain 
oscillations with an increasing amplitude. You should use the slow varying 
frequency of the perturbation and the effect of autoresonance. 

Not long ago it was found that this simple receipt is not sufficient for 
the capture into autoresonance. It was numerically obtained that there is a 
threshold value of the perturbation when the capture into resonance is not 
observed for subharmonic modes [12], [T3]. Later the similar results were 
obtained for an oscillator with a quadratic form with respect to r for the 
phase of the perturbation [13], [15]. 

In this article we have obtain a new result for two-dimensional system 
of primary resonance equations. It was shown that the autoresonant phe- 
nomenon appears when the amplitude of the perturbations is greater than 
the threshold value. This threshold value of the perturbation was found 
explicitly. 

The mathematical statement of the problem on the existence of threshold 
is reduced to an existence of asymptotic solutions in the form of power series 
with respect to as t — > 00. The general approach to the construction 
of power asymptotics can be read in [T7], [T5]. However the direct using 
of those methods does not allow one to predict the existence of threshold 
phenomenon in autoresonance equations. 



2 



1 Statement of the problem and result 

The long time evolution of amplitudes of oscillating solutions of (pQ) is reduced 
to the system of primary resonance equations: 



Our goal is to study the behaviour of the solutions of §B) when t ap- 
proaches infinity. 

It is shown there are increasing and bounded solutions. The bounded 
solutions are studied in detail. It has been determined that periodic pertur- 
bation of a system of oscillators leads to the capture into resonance. The 
asymptotic description and numerical simulations of the phenomenon are 
presented. An explicit formula for the threshold value of the perturbation 
has been found. It was found there exist solutions related to autoresonance 
phenomenon when |/| > 12. 

2 Asymptotic reduction to the system of pri- 
mary resonance equations 

In this section we present a formal asymptotic reduction of a system of pri- 
mary resonance equations. Let us construct the solution for ([1]) in the com- 
plex form 




(2) 



x = A{t) exp{iu8} + c.c, 



y = B{t) exp{2iu9} + c.c. 



(3) 



Substitute © into (P 



It yields 



A!{r) 




expjiar 2 } + ieA 



B\t) 



Neglecting the terms of order e and substituting 



A = a(r) exp{iar 2 } + c.c, B = b(r) exp{2iar 2 } + c.c. 



(4) 
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we obtain 



a (T) = —Ziara ab , 



b'(r) 



-Marb — % —^a 2 . 



Let us simplify system (jSJ). Substitution 

a{r) = XA{t), b(r) = KB(t), r = X t, 

where 



uj-Jot / a 
K — , A = uj 



Oil 



, X = - 



gives system (j2J) with / = — 1 ^ . 



(5) 



3 Numerical simulations 

In this section we present results of numerical simulations for system (T5]). 
Figures 1 and 2 show \A(t) \ for the captured and noncaptured solution. Initial 
data were identical for both cases 

4(100) = 102.669 - ^793.88, 5(100) = 386.825 + zl01.831. 

We have used / = 11.9 for the first figure and / = 12.1 for the second one. 
The second figure corresponds to the captured solution where the amplitude 
of the perturbation is greater than the threshold value. 
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Fig.l: noncaptured solution Fig. 2: captured solution 

Figures for |-B(t)| look similar for both cases. The presented numerical 
results show that there exists the threshold value of amplitude of perturbation 
that leads to essential change of structure of solution. 
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4 Algebraic asymptotic solutions 

In this section we construct algebraic asymptotic solutions of system (J2J) in 
the form 

oo oo 

A{t) = ]T a k r\ B(t) = £ b k r k , t - oo. (6) 
fc=-i fc=-i 

Theorem 1 When t —>■ oo t/iere exists the solution of system (TJ|) o/ i/ie /orm 

^(t) = -5*- 1 + ^r»+(^-^)r 6 + 0(O, 

When \ f\ > 12 and £ — > oo £/iere exzsi solutions of (EJ] o/ £/ie /orm 



= -8^cos(#)+isin(#)J* + ~r 1 + 0(* 

= -4(cos(2^) + 2sin(2^))t + (-j-2z)t" 1 + 0(r 

12 
7" 



3> 



19 

/iere sin(^) 

A 3 (t) = 8(cos(^) + isin(^))t + ^r 1 + 0(t 

f 24 

S 3 (£) = -4(cos(2^)+isin(2^))t+ (-cos(^)[^- + — ] +22[l + sin 2 (^)])r 1 



+ 0(t 



-3\ 



/ J \ 12 

nere sm(WJ = — — . 

The proof of the Theorem consists in the construction of the asymptotic 
solution of (|2J) and using Kuznetsov's Theorem |19j . 

Let us substitute representation (jSJ) into ([2]) and gather the terms of the 
same order of t. It gives the recurrent system of equations for coefficients 
<2fc, b k of asymptotic expansion (J6]) 

2ia fe + - (a^6_i + alj&fc) = (fc - 2)a fe _ 2 ~ ^ a m^> 

m,/ 

4z6 fc + -a fc a_i = (k - 2)b k _ 2 - a m a;, (7) 



here m + I = 1 — k and m,l k. 
Terms under t 2 give 



2a_i + -a*_ 1 b-.\ 



0, 46_i + -aii 



0. 



Coefficient a_i is the solution of equation 

1 



32 



a_i a_i — 2a_i = 0. 



64 and one bounded 
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It shows there are two increasing solutions with |a_i 
solution with a_i = 0. The leading-order term 6_i - 

Here we give the procedure of construction of increasing solution that 
goes to +oo when t also increases. Another increasing solution that goes to 
— oo will be discussed below. 

We suppose that 

a_i = 8exp{i^}, &_! = -4exp{2z^}. 

Terms of order of t give a homogeneous system for and &o- The matrix 
of the system for real and imaginary parts of a and b has the form 



4cos(\I>) sin(^) -4cos(\l/) 2 4sin(\I/) -4cos(*) 

4sin(\l/) 2 -4cos(\&) sin(^) 4cos(\I>) 4sin(*) 

-4sin(\I/) -4cos(*) -4 

4cos(\I>) -4sin(\&) 4 



(8) 



The rang of the matrix equals three and solution depends on a parameter. 
The solution is 



Y = HqYq = / u (sin(^), -cos(^), -sin(2^),cos(2^)). 

It yields 

a = fi \ sin(^) - % cos(^)| , b = [i \ - sin(2\I>) + i cos(2^) 



Terms of order of t° give a nonhomogeneous system for real and imaginary 
parts of ai and b\ with the matrix (jSJ). To obtain a non-trivial solution the 
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right-hand side 

F = ( - 8 cos(#) - X -nl sin(tt), -/ + ^ cos(#) - 8 sin(tt), 

4 cos(^) 2 - cos(^) sin(^) - 4 sin(^) 2 , 

cos(^) 2 + 8 cos(^) sin(^) - i/i 2 , sin(^) 2 ) 

should be orthogonal to the solutions of the adjoint system. The solution of 
adjoint system is 

Z = (-cos(tf),-sin(tf),cos(2tf),sin(2tf))- (9) 

The solvability condition for the system for real and imaginary parts of a\ 
and b\ has the form 

sin(tf) = -j. (10) 

The variable \1/ that determines a turning of the leading-order term of asymp- 
totic expansion can be determined when |/| > 12. 

Real and imaginary parts of cii and b\ are represented as a sum of general 
solution for homogeneous equations and partial solution of nonhomogeneous 
equation. 

(K[oi],9f[oi],3fi[6i],9f[6i]) = + 0, 

cos(fr)(-192 - 2 J 2 + j> 2 cos(tt)) 

8/ ' 
1152/ cos(*) + 8/ 3 cos(^) + 864^ 2 - 9/V 2 xT 



4/ 3 cos(^) 



') 



Coefficients for higher-order terms of asymptotics are determined in the sim- 
ilar way. The solutions for these coefficients depend on a parameter This 
parameter is determined from the solvability condition for the nonhomoge- 
neous system of algebraic equations with degenerate matrix. The solvability 
condition looks as follows 

Z-(F 1 ,F 2 ,F 3 ,F 4 ) = 0. (11) 

where Z is a solution of the adjoint system, (iq, F 2 , F3, F4) is the right-hand 
side of equation for real and imaginary part of ak+2, bk+2- 
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Vector (F 1} F 2 , F 3 , F4) contains a product of coefficients for higher-order 
correction terms of The algebraic system is written for 



ft[a fc ],9[a fc ], &[&*], 3[6 fc ] 

It yields the following rule for right-hand side 

F 1 = -3?[a fc _ 1 ] + 3?[^a^ 
F 2 = -9f[a fc _i] + %[^2a* m bi 
F 3 = -&[b k - 1 ]+&[£ l a m a l 

F A = -9f[6 fc _i] + 3 [ a mCll 



(12) 



The parameter fi k is determined from relation (TTTTl for all values of fc. The 
equation for \i k looks as follows 



z-(y ©x fc )-£-x, 



0, 



(13) 



where Z is a solution of homogeneous union system, Y is a solution of ho- 
mogeneous system and X k is a partial solution of nonhomogeneous system 
for k-th correction term. An operation is determined by 



YQX 



+ va%2 + yix 3 + 1/2x4 
-y-i'x 2 + y 4 xi - y A x z + 2/ix 4 

2j/ia;i - 2y 2 x 2 
2y 2 xi + 2j/ix 2 



We determine the parameter \x k and relation fflBl becomes valid. 
The solution X k can be expanded on a basis Y , Yi, Y 2 , Y 3 , where 

Y = (sin(^), -cos(^), -sin(2^),cos(2#)V 

Y x = (cos(^),sin(*),0,o), 

Y 2 = (o,0,cos(2tf),sin(2tf)), 

y 3 = (o,0,-sin(2^)cos(2^),cos(2^)sin(2^)). 

Direct calculations give 

Z ■ (Y Yi) = 0, for z = 0,l and Z ■ (Y Y) ^ 0, for i = 2, 3, 



(14) 
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and the vector Z is not orthogonal to vectors Yi,Y 2 ,Y 3 . It yields a nontrivial 
equation for \i k 

fi k z ■ (y © \c K2 y 2 + c K3 y 3 }) = z ■ (c k . 1>x Y x + c fe _ 1;2 r 2 + c fe _ 1)3 y 3 

Note. We can construct another solution with the leading-order term 
a\ = — 8 in similar way. The change of sign for a\ does not lead to essential 
change of procedure of asymptotic construction. 

Above asymptotic constructions give the following algebraic asymptotic 
expansions for increasing solutions 

A 1 (t) = -8(cos(^) + zsin(*))t + ^r 1 + 0(r 3 ), 

£i(t) = -4(cos(2^) + 2sin(2^))t + (-^-2z)r 1 + 0(r 3 ), (15) 

12 

where sin(^ r ) = — . 

A 3 (t) = 8(cos(^) + isin(^))t+ ^r 1 + 0(r 3 ), (16) 

f 24 

B 3 (t) = -4(cos(2^)+isin(2^))t+ (-cos(^)[^ + y] +2z[l + sin 2 (^)])^ 1 

+ 0(r 3 ), (17) 

12 

where sin(^) = — — . 

We construct the finite solution of the form 

oo 

A(t) = x> fc r fc , 

k=l 

oo 

B(t) = E^. (18) 

k=l 

Substitution of (|T8|) into equation gives the recurrent sequence of problems 
(JTj) for coefficients of ({TBI . The coefficients a k and b k are determined from a 
system of algebraic equations with a nondegenerate matrix. 

Relations of order of t° give the following system for real and imaginary 
parts of the leading-order terms 



2ia± r — 2au = —if, 4ib\ r — 4b\ 



i=0- 
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These equations allow us to determine the leading-order term of asymptotic 
expansion 

/ 

ai = --, bx = 0. 

The real and imaginary parts of 0,2, &2 should equal zero. 
Relations of order of t~ 2 give nontrivial equations 

/ if 2 
2ia 3r - 2a 3i = --, Aib 3r - Ab 3i = - — . 

2 lb 

It yields 

= if h . = Jl 
Qs 4 ' 3 64 ' 

The recurrent procedure gives 

Mt) = -^r 1 + ^r 3 + (^-^)t- 5 + 0(r 7 ), 

f 2 7? f 2 

s 2 m = -^-r 3 + ^-r 5 + o(r 7 ), (19) 

v ; 64 256 K J K ' 



5 Neighborhoods of equilibrium positions 
5.1 Stability by linear approximation 

In this section we present the analysis of stability for the constructed alge- 
braic asymptotic solutions by linear approximation. 

Let us consider the system of equations that are linearized on the increas- 
ing asymptotic solution (Ai(t) , Bi(t)) . Substitute 

a(t) = Ai{t) + a(t), bit) = B x {t) + p{t), 

into (j2J). It gives the system of linear differential equations for a(t) and (3(t). 
It is convenient to rewrite this system as a system for real and imaginary 
parts of a(t) and 0(t). The eigen- values for the system are 

\ 1 = -Aiy/St + Oir 1 ), A 2 = 4iv / 3t + 0(r 1 ), 

4/ p _ 144 x/f 1 - 144 

a 3 = - V/ ^ + or 1 ), a 4 = v/ ^_ + or 1 ). 
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A4 has a positive real part when / > 12. It yields that the constructed above 
algebraic solution is not stable with respect to small perturbations. 

The similar linearization on the algebraic solution A2(t),B2(t) leads to 
the matrix with eigen- values of form 

\ 1 = -ut + o(r l ), \ 2 = 4it + 0(r 1 ), 

A 3 = -2it + 0(r l ), \ A = 2it + 0(r 1 ). 

The leading-order terms of the given asymptotic expansions are imaginary 
and a question on stability of A 2 (t),B 2 (t) requires additional investigations. 

A matrix of the system of linear differential equations that are realized 
on A3 (t) , B$ (t) has eigen- values of the form 

\ 1 = -Ut + 0(t _1 ), \ 2 = Ait + 0(r 1 ), 



A3 = -* ^= + 0(t ), A 4 = z + 0(t ). 

It leads to additional investigations of stability of A 3 (t),B 3 (t) for / > 12. 



5.2 Oscillating asymptotic solution in the neighbor- 
hood of bounded solution 

In this section we present an oscillating asymptotic solution in the neighbor- 
hood of bounded solution. The constructed solution of (j2J) has the form 

A(t) = a{t) exp{-it 2 }, B{t) = b{t) exp{-2it 2 }. (20) 

The substitution gives 

ia'{t) = -a*b + f exp{it 2 }, 

ib'(t) = \a\ (21) 

Here we study the behavior of solution of the system in the neighborhood of 
bounded asymptotic solution (A 2 ,B 2 ). Substitution 

a = A 2 exp{it 2 } + a, b = B 2 exp{2zt 2 } + (3, (22) 
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$ives a system for a, (3 

1 1 



id = ^ a *(3 + ^(A* 2 f3 exp{-it 2 } + a* B 2 exp{2it 2 }), 



i[3' = I« 2 + ^A 2 aexp{it 2 }. (23) 



Theorem 2 There exists a formal asymptotic solution of (TJ|) in the form of 
{2E with 



a = J> fc r fc , (3 = J2Pkt- k . (24) 

fc=0 k=0 

This asymptotic solution depends on four real parameters. The leading-order 
terms of the expansion are determined in terms of elliptic functions. 

The same result for the system of three coupled oscillators was obtained in 

The proof of the theorem will be obtained by formal asymptotic con- 
struction of flU]). Substitute (J22J, (ED into (EH} and take into account 1115]) . 
Gathering the terms of the same order of r gives a recurrent system of equa- 
tions for coefficients of the asymptotic expansion. The leading-order terms 
are satisfied to 

id (t) = ^a* /3 , 

i(3' (t) = \a 2 . (25) 

This system can be solved in terms of elliptic functions. The system has two 
conservation laws 

\a \ 2 + 2\p \ 2 = E 2 , 
K) 2 /? + (a ) 2 {3* = H. 

iH 

Function is a Hamiltonian of (1231) . The first conservation law allows us 

4 V! y 

to obtain the following representation for ao> A) 

E 

«o = Eexp{iip} cos(\l/), /3q = — j= exp{iijj} sin( 1 J r ). 

V2 
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After this substitution the second conservation law becomes 

H = v / 2 J B 3 cos 2 (^)sin($)cos($), where $ = 2y? - ip. (26) 
System ( 1251) can be written as 

V' = ~^|COs($)siii(tf), 

V>' = r cos($)cos 2 (^)sin" 1 (^), (27) 

2y 2 

*' = 2V2 sin ( $ ) C0S W' 

Using expression fl26l) allows us to separate system (!27j) . It yields the separate 
equation for \& 



v 2£ 6 cos 4 (#)sin 2 (^) - H 2 
4E 2 cos(^) sin(^) 

The solution of the equation is 



/■cos(2<I') W 7 i 

/ ,r f f= = -T' (28) 

Jmo v G — ir — ir + it l 
E 6 - AH 4 



where G = — . Function cos (2$) is a bounded periodic function with 

respect to t. We determine the function $ from (|2"6"|) and integrate the first 
equation of (J27J) 

H r* dt 



AE 2 J l + cos 2 (2^) 

The integrand is a periodic function with respect to t with nonzero average 
value. It yields the following behaviour ip — 0(t), t —>■ oo. The function ip is 
determined by ip = 2<p — $. Thus we have constructed a family of solutions 
chq, (3q that depends on four parameters E,G,u ,tp . 

The first correction terms of expansion f[2"2"j) are determined from the 
linearized system 

1 f 

ia[ - -{alfio + otlPx) = -jPo exp{-it 2 }, 

ij3[ — -aocti = — — «o exp{it 2 }. (29) 
13 



The fundamental matrix W of the homogeneous system that relates to (I2"9j) is 
formed by the first order derivatives with respect to parameters of solutions 
for fl25l) . When t —>■ oo we obtain the following asymptotic behavior 

d E a = 0(t), d G a = O(t), d U0 a = O(l), d^ao = O(l). 

The similar formulas are valid for functions ckq, (3q, /3q- Thus the fundamental 
matrix contains two columns of the order of t and has the determinant of the 
order of a constant. 

The solution of the nonhomogeneous system is increasing due to fast 
oscillations of the right-hand sides. It is convenient to rewrite the right-hand 
side of the system as a sum of two vectors g + exp{zt 2 } + g~ exp{— it 2 }. Here 



0) 



fPo fao 



4 4 

After integrating twice we obtain 



0, 
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(-)^00 



w 



+ Tr 

+ 



W' 1 g + exp{it 2 }dt = |^-exp{it 2 } 

■W~ 



W.W^g- 



t 



g 



it exp{ ^ } + 1 

W~ 1 g+ 



pod r/ i/i/ ~ ± n~ r \ i \ -i i 
'exp{*t 2 } - Wj t [(— 2/-) ^] exp{^^ 2 }d^ 



t )'exp{it 2 }+ f 



2t 
W{W~ 



1 )'g + 



exp{it } 



W 



2t 



* 1 

2t 



exp{it 2 }dt. 



RangW = 2 and any minor of the matrix does not contain terms of the order 
of t 2 . Taking into account that the determinant of matrix W is of order of 
constant we obtain the order of terms of inverse matrix W~ l . Their order is 



t. Thus term 



has the order of a constant. Another integral can 
t t _ 6 

be evaluated by parts and estimated by 0(t 1 ). The similar evaluations are 

valid for the part of the solution with g~. It yields that the first correction 

terms a±, Pi are finite. 

The next order correction terms are determined from the system of type 

of (1291) . The right-hand sides of the systems are quadratic forms of previous 

corrections /3k and coefficients of (IT31) . Solutions are bounded. The 

theorem is proved. 
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5.3 A neighborhood of the increasing algebraic solu- 
tion 

We investigate a neighborhood of the increasing asymptotic algebraic solution 
(A 3 (t) , B 3 (t)) numerically. We solve an initial value problem for big values 
of t. Initial data are close to the algebraic asymptotic solution. Denote by 
A 3 (t;2) and B 3 (t;2) a sum of two first terms of the asymptotic expansion 
(A 3 (t),B 3 (t)). We choose / = 12.1 and 

A\ t =i 00 = A 3 (100; 2) + 0.1, S| t=100 = #3(100; 2) + 0.1; 

and solve system ([2]) numerically. 

The results are presented below. The parameter t goes from 100 to 150. 
Here we graph on the complex plane A and B. 




Fig. 3: A 3 (t) is a straight line and Fig.4: B 3 (t) is a straight line and 
numerical solution A(t) numerical solution B(t) 

Here we present the graph for comparative difference for numerical solu- 
tion and algebraic asymptotic solution. 
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Fig.5: comparative 
difference: 



A(t)-A 3 (t) 
A(t) 



110 120 130 140 150 



Fig. 6: Absolute value of 
numerical solution 



The graph of numerical solution \B(t)\ looks the same. The constructed 
above numerical solution proves the existence of oscillating solutions near the 
increasing algebraic asymptotic solution. 

We thank Prof. L.A. Kalyakin for stimulating discussions and Prof. A.D. 
Bruno for discussions. We are also grateful to Anna A. Shchipitsyna, PhD 
for helping us to prepare the presentation of our results. 
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